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NUMERICAL  HYDRODYNAMIC  CALCULATIONS  OF  THE  FLOW  OF  THE 
DETONATION  PRODUCTS  FROM  A  POINT-INITIATED  EXPLOSIVE  CYLINDER 

By 

D.  Piacesi,  Jr. 

ABSTRACT:  With  the  use  of  the  NOL  two-dimensional  hydrodynamics 
computer  program  (CYCLONE) ,  a  detailed  numerical  analysis  is  made 
of  the  flow  in  the  product-gas  of  a  detonating  cylinder  of  pentolite 
(50/50  PETN/TNT)  explosive.  The  explosive  cylinder  is  initiated  at 
a  point  on  the  central  charge  axis  at  one  end- face.  The  detonation 
front  is  assumed  to  be  spherical.  This  is  consistent  with  experimen¬ 
tally  determined  wave  shapes.  The  calculations  are  carried  out  to 
the  point  where  the  detonation  front  is  at  1.71-charge-diameters  from 
the  point  of  initiation.  The  results  of  the  calculations  show  that 
the  flow  in  the  "detonation  head"  (region  bounded  by  the  detonation 
front,  and  the  lateral  and  rear  rarefaction  waves)  is  entirely 
spherical.  Moreover,  this  spherical  flow  is  exactly  that  which  is 
described  by  the  similarity  solution  of  G.  I.  Taylor  for  a  spherical 
detonation  wave. 
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NUMERICAL  HYDRODYNAMIC  CALCULATIONS  OF  THE  FLOW  OF  THE 
DETONATION  PRODUCTS  FROM  A  POINT-INITIATED  EXPLOSIVE  CYLINDER 

This  report  discusses  the  results  of  numerical  hydrodynamic 
calculations  on  the  flow  in  the  product-gas  resulting  from  the 
detonation  of  a  cylinder  of  explosive.  This  report  will  be  of 
interest  to  scientists  who  are  dealing  with  experiments  involving 
point-initiated  explosive  cylinders,  and  to  engineers  concerned 
with  conventional  explosive  weapons  design.  This  work  was  carried 
out  under  NAVORD  Task  ORD-033-221/092-1/F008-08-H  Problem  004  - 
Explosion  Hydrodynamic  Calculations. 
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I.  INTRODUCTION 

Let  us  examine  the  flow  produced  by  the  detonation  of  an 
unconfined  right  circular  cylinder  of  high  explosive,  where  the 
initiation  is  at  a  point  on  the  central  axis  at  one  end-face 
(refer  to  Fig.  1) .  As  the  detonation  progresses  from  the  point  of 
initiation,  the  gaseous  explosion  products  expand  out  through  the 
end-face.  Consequently,  a  rarefaction  wave  travels  into  the  gas 
behind  the  detonation  front.  When  the  detonation  reaches  the  side 
of  the  cylinder,  a  shock  is  transmitted  into  the  surrounding  medium 
and  a  lateral  rarefaction  wave  is  reflected  into  the  explosion  products 
Subsequently,  the  head  of  the  lateral  rarefaction  wave  impinges  on  the 
central  charge  axis  and  annihilates  the  rear  rarefaction  wave.  The 
region  bounded  by  the  head  of  the  lateral  rarefaction  and  the 
detonation  front  is  sometimes  referred  to  as  the  "detonation  head".1* 

While  considerable  attention  has  been  given  to  the  detonation 
head,  there  does  not  appear  to  be  in  the  literature  any  attempt  to 
calculate  the  actual  flow  in  the  detonation  head  region.  With  the 
use  of  the  Naval  Ordnance  Laboi  n  U. ..  y  '  a  two-dimensional  hydrodynamics 
computer  program,  CYCLONE^ ,  a  detailed  numerical  analysis  of  this 
transient  compressible  flow  problem  can  be  carried  out.  The  results 
of  calculations  for  a  two-inch  diameter  cylinder  of  pentolite  (50/50 
PETN/TNT)  explosive  are  presented  in  Sec.  V  of  this  report;  The 
calculations  were  performed  on  an  IBM-7090  computer. 

*  References  are  on  Page  15. 
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II.  THE  PHYSICAL  MODEL 

We  make  the  physical  assumptions  about  the  detonation  that 
(1)  the  chemical  reaction  zone  is  infinitely  thin  (i.  e.,  the  solid 
explosive  is  instantaneously  converted  to  gaseous  products  as  the 
shock-detonation  front  passes  over  it) ,  (2)  the  Chapman- Jouguet  (CJ) 
conditions3  hold  at  the  detonation  front,  (3)  the  detonation  propagates 
in  all  directions  from  the  point  of  initiation  with  the  CJ  detonation 
velocity,*  (4)  the  explosive  is  surrounded  by  a  vacuum,  and  (5) 

"edge"  effects  are  negligible.  Here,  the  term  "edge"  effects  refers 
to  the  radical  change  in  the  curvature  of  the  detonation  front  near 
the  cylindrical  surface  of  the  charge.  This  effect  is  produced  by 
the  influence  of  the  lateral  rarefaction  wave  on  the  chemical  reaction 
zone. 

In  solid  high  explosives,  the  width  of  the  chemical  reaction  zone 
is  of  the  order  of  several  millimeters  and  therefore  assumption  (1) 
does  not  impose  a  too  unreal  restriction  on  the  physical  model. 
Assumption  (1)  does,  however,  eliminate  any  effect  of  the  lateral 
rarefaction  wave  on  the  chemical  reaction  zone  at  the  sides  of  the 
cylinder.  Hence,  assumption  (5)  results  as  a  natural  consequence  of 
(1).  In  the  experiments  by  Cook,  et  al.1,  the  edge  effect  did  not 
exceed  2  mm  in  any  case.  Assumption  (2)  is  introduced  into  the 
calculations  via  the  equation  of  state  of  the  detonation-product  gas 
( see  Sec .  IV) . 

_ _ _  i 

*  The  CJ  detonation  velocity  is  defined  as  D  =  \^  L  (pa  -pi  )/(\*  )  ]5' 

where  pi  ,  Vi  (p  --  pressure,  v  «=  specific  volume)  is  the  initial  state 
of  the  solid  unreacted  explosive  and  ps  ,  V2  is  the  point  of  tangency 
of  the  Rayleigh  line  from  Pi  ,  Vi  to  the  reactive  Hugoniot  of  the  product 
gas.  This  is  the  theoretical  maximum  steady-state  velocity  which  the 
detonation  can  achieve.  See  Ref.  (3). 
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The  shape  of  the  detonation  front  directly  affects  the  flow  in 

the  region  behind  the  wave  front.  Assumption  (3)  imposes  the 

requirement  that  the  wave  front  coincide  with  the  surface  of  an 

expanding  sphere,  whose  radius  of  curvature  is  equal  to  the  product 

of  the  CJ  detonation  velocity  (DOT)  and  the  time  (t)  elapsed  after 

initiation.  Cook,  et  al.,1  and  Jaffe  and  Clairmont4,  show  experimentally 

that  in  point-initiated  cylinders  of  "ideal"  explosives*  the  detonation 

front  is  indeed  initially  spherical,  and  that  the  radius  of  curvature 

of  the  front  is  D__t.  However,  there  is  a  fundamental  difference 

CJ 

between  the  two  sets  of  data.  In  the  work  by  Cook  et  al.,  the 
radius  of  curvature  of  the  detonation  front  approaches  a  constant 
value  by  the  time  the  wave  has  propagated  a  distance  of  less  than 
four  charge  diameters  from  the  point  of  initiation.  From  this  point 
on,  the  detonation  front  propagates  in  a  "steady  state"  condition. 

In  the  experimental  work  by  Jaffe  and  Clairmont  with  two-inch  diameter 
tetryl  cylinders,  the  radius  of  curvature  of  the  detonation  front 
increases  geometrically  (radius  =  D  t)  at  a  distance  of  f^ur  charge 
diameters  from  the  point  of  initiation  with  no  indication  of  reaching 
a  steady  state.  In  the  computer  calculations  which  are  described  in 
this  report,  the  mathematical  model  of  the  detonation  is  consistent 
with  both  sets  of  data,  since  the  calculations  were  only  carried  out 
until  the  detonation  front  is  at  1.71  charge  diameters  from  the  point 
of  initiation.  There  is,  however,  an  obvious  extrapolation  of  the 
calculations  which  can  be  made  in  the  light  of  the  ja f fe-Clairmont 
experimental  results.  This  is  done  in  Sec.  VI. 

*  We  define  an  "ideal"  explosive  as  one  in  which  the  detonation 
propagates  at  the  theoretical  maximum  steady-state  velocity  (i.e.,  the 
CJ  detonation  velocity). 
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III.  THE  HYDRODYNAMIC  EQUATIONS  AND  THE  MATHEMATICAL  MODEL 
We  can  write  the  hydrodynamic  equations  for  transient,  isentropic, 
axially  symmetric  flow  in  Lagrangian  form  as  follows:  Denote  the 
Lagrangian  coordinates  by  k,  1,  the  Eulerian  coordinates  by  R,Z, 
the  velocity  components  by  u,  v,  the  density  by  n,  the  specific  volume 
by  V  =  (1/p),  the  internal  energy  by  E,  the  pressure  by  p,  and  the 
time  by  t.  Then  the  equations  are 


(continuity) 


where 


(momentum) 


R 


RR 

R„R0  ' 

3(M)  ' 


d‘R  -  -R 

3E 

PZ 

dt"  n  R/to 

3  k.  U 

Pi  3k 

dJ2  -  -R 

dtl  1 

Pi  Pk 

Pi.  H_ 

) 


) 


(1) 


(2) 


(4) 


(energy) 


7E  _  u  3V 


(5) 


where  the  subscript  zero  refers  to  the  initial  state  of  the  material. 
The  equation  of  state  for  each  material  is  assumed  given  in  the  form 


p  =  p  (E ,  V)  . 

The  von  Neumann-Richtmyer  finite  difference  method,'  which 
handles  shocks  automatically,  is  used  for  the  numerical  integration 
of  the  hydrodynamic  equations.  In  this  scheme,  an  artificial 
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dissipative  term  is  introduced  into  the  conservation  of  energy  and 
momentum  equations.  This  term  causes  shocks  to  have  a  finite  thick¬ 
ness  and  makes  the  flow  variables  change  abruptly  but  continuously 
through  the  region  of  the  shock.  In  addition,  the  correct  entropy 
change  is  effected  across  the  shock.  In  CYCLONE,  the  artificial 
dissipative  term  is  defined  as 


a  -  j>  A  /W  f 
^  v  v/  Ut  J 


(6) 


where  the  variables  A , V ,  and  t  are  defined  above  and  b  is  the  shock- 
width  constant,  which  can  be  adjusted  to  give  the  desired  sharpness  to 
the  shock.  The  pressure,  p,  in  equations  (3),  (4),  and  (5)  is  then 

replaced  by  P  =  (p+q).  The  finite  difference  analogues  of  the 
hydrodynamic  equations  which  are  actually  used  in  the  numerical 
integration,  can  be  found  in  Ref.  2. 

Let  us  now  briefly  examine  the  mathematical  model  that  is  used 
in  the  finite  difference  method.  As  a  result  of  the  cylindrical 
symmetry,  the  three-dimensional  flow  can  be  represented  by  the  flow 
in  any  R,  Z  (R^O,  Q  =  const)  half  plane  (Fig.  2).  Lagrangian 
coordinate  (k,i)  lines,  imbedded  in  the  material,  divide  the  material 
into  zones.  Initially,  the  mass  of  each  zone  is  calculated.  This 
mass  remains  constant  throughout  the  computations,  thus  assuring 
conservation  of  mass.  There  is  then  assigned  to  each  mesh  point 
(intersection  of  a  k  and  an  1  line)  a  mass  equal  to  one-fourth  of  the 
mass  in  the  tour  surrounding  zones.  The  velocity  components  of  each 
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mass  point  and  the  pressure,  internal  energy,  and  density  of  each  zone 
are  given  as  initial  conditions. 

At  each  computation  cycle,  the  program  calculates  consecutively 

for  each  zone  p  from  the  continuity  equation,  E  and  p  from  the 

simultaneous  solution  of  the  energy  equation  and  the  equation  of  state 

and  q  from  the  formula  given  in  eq.  (6).  At  this  point,  a  time  step 

0 

for  the  integration  is  determined  from  stability  conditions.  The 
accelerations  of  each  of  the  mass  points  are  next  determined  by  the 
momentum  equations.  From  the  accelerations  and  the  time  step,  the 
new  velocity  components  and  positions  of  the  mass  points  are  then 
obtained.  The  calculations  are  advanced  in  finite  increments  of  time 
by  repeating  the  above  steps  at  every  integration  cycle.  The  material 
flow  is  thus  determined  by  the  motion  of  the  mass  points. 

A  programmed  detonation  scheme  is  used  to  instantaneously 
convert  the  solid  high  explosive  to  gaseous  explosion  products.  Since 
we  assume  the  detonation  propagates  with  the  constant  CJ  velocity 
from  the  point  of  initiation,  we  therefore  can  calculate  the  position 
of  the  detonation  surface  at  each  computation  cycle.  For  each  zone 
containing  solid  explosive,  a  test  is  made  to  see  if  the  detonation 
front  has  arrived  at  the  zone  (corner  nearest  the  point  of  initiation) 
If  it  has  arrived,  the  chemical  energy  of  the  solid  explosive  which 
is  contained  in  the  zone  is  added  to  the  internal  energy  of  that  zone* 
For  this  and  subsequent  time  cycles,  this  zone  uses  the  equation  of 

state  for  the  gaseous  explosion  product:-:.  If  the  detonation  has  not 

*  In  the  reported  calculations  (Sec.  V),  the  chemical  energy  was 
actually  added  in  discrete  amounts  of  1/3  the  total  chemical  energy 
in  the  zone  as  the  detonation  front  passed  over  each  1/3  section  of 
the  diagonal  of  the  zone.  By  allowing  the  energy  adjustment  to  take 
place  over  several  computation  cycle-.,  the  typical  oscillations  found 
in  numerical  calculation-',  of  this  type  were  greatly  reduced. 

6 


NOLTR  66-150 


yet  arrived,  the  equation  of  state  of  the  solid  unreacted  explosive 
is  used. 

IV.  EQUATIONS  OF  STATE 

In  compressible  fluid  flow  calculations,  it  is  necessary  to 
have  an  equation  of  state  for  each  material.  The  equations  of  -tate 
which  are  used  for  the  solid  unreacted  explosive  and  the  detonation 
product  gases  are  discussed  in  this  section. 

When  making  numerical  hydrodynamic  calculations,  it  is  often 
a  satisfactory  approximation,  at  least  for  solids,  to  use  the 
shock  Hugoniot  (the  locus  of  possible  states  which  are  attainable  from 
some  initial  state  by  a  single  shock  transition)  as  an  equation  of 
state.  The  equation  of  state  for  the  solid  unreacted  Pentolite 
explosive  is  taken  as 

p  <*  0.1282  it  +  1.193  n3  ,  (7) 

where  p  is  the  pressure  (megabars)  ,  u  =  p/p  -1,  and  p  is  the 

O 

density  (g/cm3).  The  subscript  zero  refers  to  the  initial  state. 
Equation  (7)  may  not  closely  represent  the  shock  Hugoniot  for 
unreacted  pentolite.  It  is  probably  more  representative  of  the 
shock  Hugoniot  for  unreacted  Composition  B  explosive  (see,  for 
example,  Ref.  (7)).  However,  almost  any  reasonable  form  for  the 
equation  of  state  is  adequate  for  our  calculations  since  the 
detonation  is  programmed  as  a  function  of  time  and  does  not  depend 
on  the  state  of  the  solid  explosive.  The  form  of  the  equation  of 
state  for  the  solid  explosive  has  essentially  no  effect  on  the  flow 
in  the  detonation-product  gases. 
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The  equation  of  state  of  the  pentolite  detonation- product 
gas  is  assumed  to  be  given  in  the  form8 

pa  (Ap  +  Bp2)  E  +  CP3  (8) 

where  A,  B,  and  C  are  constants  and  p,  E,  and  p  are  the  thermo¬ 
dynamic  variables  described  above.  The  value  A  =  0.35  is  obtained 
by  defining  A  =  (7T-1),  whereYis  the  average  ideal  gas  specific 

heat  ratios  (c  /c  )  of  the  detonation  products.  Equation  (8)  thus 
p  v 

approaches  the  polytropic  gas  equation  for  small  p.  The  constants 
B  =  0.1243  and  C  =  0.01279  were  chosen  so  that  eq.  (8)  would  reproduce 
the  particular  CJ  detonation  conditions  and  the  CJ  isentrope  obtained 
by  thermochemical  calculations.* 

The  CJ  detonation  conditions  for  pentolite  (p  =  1.65)  gotten 

O 

with  the  thermochemical  calculations  are  pressure  pCJ  =  0.2452 
megabars,  density  p^T  =  2.210  g/cm3 ,  sound  speed  c^^  =  0.5714 
cm/psec,  and  detonation  velocity  D  =  0.7655  cm/psec.  The  value 
of  the  CJ  specific  internal  energy  E  =  0.0775  megabar-cm  /g 
is  obtained  by  calculating  the  area  under  the  CJ  isentrope  from 
the  CJ  pressure  to  zero  pressure.  The  chemical  energy  Eq  which  is 
used  in  the  computer  calculations  is  obtained  from  the  Hugoniot 


relation 


E  -  E  =  (P  +  p  )  (V  -  V  T ) / 2 . 
CJ  0  CJ  O  O  CJ 


This  gives  E  =  0.0587  megabar-cm3 /g .  E  is  the  energy  which  is 


added  to  the  specific  internal  energy  of  a  zone  as  the  detonation 
front  reaches  the  zone  (see  Sec.  III). 

*  The  thermochemical  calculations  were  carried  out  by  H.  Hurwitz, 
using  the  RUBY  computer  program9  which  contains  the  Cowan-Fickett 
version  of  the  Becker-Kistiakowsk^-Wilson  equation  of  state  .  The 
modified  RDX  parameters  of  Maderirwere  used  in  the  calculations. 
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V.  RESULTS  OF  CALCULATIONS 

The  results  of  the  calculations  on  the  flow  produced  by  the 
detonation  of  a  two-inch  cylinder  of  pentolite  explosive  are 
summarized  in  Figs.  3-12.  Figures  3-10  and  12  were  processed  directly 
from  the  computer  calculations  (excluding  the  labels  and  dashed-in 
rarefaction  fronts)  on  a  Stromberg-Car Ison  SC-4020  Microfilm 
Recorder . 

The  cross  section  (parallel  to  the  central  axis)  of  the 
detonating  cylinder  is  shown  in  Figs.  3-6  for  different  instants 
of  time  after  initiation.  Initiation  of  the  detonation  is  at  t=0. 

The  Lagrangian  computation  grid  is  plotted  on  the  cross  section.  The 
ini' ■.•■<•'.1  grid  (Fig.  3)  consists  of  approximately  4000  mesh  points. 
Additional  grid  is  generated  as  the  detonation  progresses  down  the 
cylinder.  At  t  =  11.3841  microseconds  the  computation  grid  has  9000 
mesh  points  (the  top  half  of  the  grid  shown  in  Fig.  6) .  The  initial 
zone  size  of  the  solid  explosive  is  0.0508  x  0.05  cm.  Individual 
cones  in  the  solid  explosive  and  the  compressed  explosion  products 
cannot  be  clearly  distinguished  in  the  figures  because  of  the 
inadequate  resolution  of  the  scale  which  was  chosen  for  the  plots. 

In  Figs.  4-6,  however,  the  zone  configuration  is  discernible  in  the 
region  where  the  gas  has  undergone  some  expansion.  At  2.9753  micro¬ 
seconds  after  initiation  (Fig.  4) ,  the  detonation  has  not  yet  reached 
the  side  of  the  cylinder.  The  detonation-product  gas  expands  only  out 
the  end-face  of  the  cylinder.  After  the  detonation  reaches  the  side 
of  the  cylinder  (Fig.  5)  the  gas  expands  radially  as  well.  The  cross 
section  at  t  =  11.3841  microseconds  is  shown  in  Fig.  6.  The 
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Lagrangian  grid  line  which  was  originally  the  end-face  surface  of  the 
solid  explosive  was  dropped  from  the  calculations  at  t  =  9.9171  micro¬ 
seconds  and  therefore  does  not  appear  in  Fig.  6. 

The  pressure  distribution  in  the  product  gas  at  specific  times 
after  initiation  is  illustrated  in  Figs.  7-10.  Here,  the  positions 
of  the  isobars  in  the  cross  section  of  the  detonation  are  shown. 

The  leftmost  isobar  is  the  10  kilobar  line.  On  each  consecutive 
isobar  to  the  right,  the  pressure  increases  by  10  kilobars.  The 
detonation  front  is  at  the  CJ  pressure  of  245  kilobars.  The  head 
of  the  rear  and  lateral  rarefaction  waves  are  shown  by  dashed  lines. 
The  head  of  the  lateral  rarefaction  wave  can  be  located  on  each  isobar 
(s  60  kb)  at  the  point  where  the  isobar  departs  from  being  circular. 

On  the  central  axis,  the  head  of  the  rear  rarefaction  wave  moves  with 
the  53  kb  isobar.  In  Fig.  7,  only  the  rear  rarefaction  wave  exists. 

In  Figs.  8  and  9,  both  rarefactions  can  be  seen.  The  head  of  the 
lateral  rarefaction  impinges  on  the  central  axis  at  t  =  11.3841  micro¬ 
seconds  (Fig.  10) ,  thereby  annihilating  the  rear  rarefaction  wave. 

At  subsequent  times,  the  head  of  the  rear  rarefaction  wave  does  not 
appear . 

It  should  be  noticed  that  the  flow  is  entirely  spherical  in 
the  region  (detonation  head)  bounded  by  the  heads  of  the  rear  and 
lateral  rarefaction  waves  and  the  detonation  surface.  Moreover,  this 
spherical  flow  is  exactly  that  which  is  described  by  the  similarity 
solution  of  G.  I.  Taylor  for  a  spherical  detonation1^  (Hereafter 
referred  to  as  the  spherical  laylor  wave).  In  Fig.  11,  the  pressure 
and  particle  velocity  profiles  along  the  central  axis  in  the  CYCLONE 
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calculations  (t  =  7.6359  microseconds)  are  compared  to  the  pressure 
and  particle  velocity  distribution  in  a  spherical  Taylor  wave.* 

The  two  sets  of  data  are  shown  to  be  in  exact  agreement  in  the  region 
where  the  Taylor  wave  particle  velocity  is  greater  than  zero.  The 
departure  of  the  two  data  in  the  Taylor  wave  region  of  zero  particle 
velocity  is  due  to  the  two  different  boundary  conditions  at  the 
initiation  point.  In  the  Taylor  similarity  solution  the  detonation 
is  spherical  and  the  geometry  imposes  the  boundary  condition  that 
the  gas  at  the  point  of  initiation  remain  fixed.  In  the  CYCLONE 
calculations  tne  detonation  is  hemispherical  and  v.he  gas  at  the 
point  of  initiation  on  the  central  axis  is  allowed  to  expand  freely. 

The  expansion  of  the  gas  along  the  axis  causes  a  rarefaction  wave 
to  follow  the  detonation  front.  This  rear  rarefaction  wave  exists  only 
until  it  is  annihilated  by  the  lateral  rarefaction  wave  impinging  on 
the  central  axis.  Previous  to  this  time,  the  penetration  of  the 
rear  rarefaction  wave,  along  the  axis,  is  known  since  the  position  of 
the  head  of  this  wave  must  correspond  with  the  "sonic  point"  (point  of 
zero  particle  velocity)  in  the  Taylor  wave.  We  see  from  Fig.  11  that 
this  is  sc . 


*  This  Taylor  wave  was  numerically  generated  with  the  use  of 
W.  A.  Walker's  computer  routines1*1  .  The  equation  of  state  for 
the  detonation  products  given  in  Sec.  IV  was  used  in  the  Taylor 
wave  calculations. 
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The  effect  of  the  von  Neumann-Richtmyer  finite  difference  method 
on  the  flow  variables  at  the  shock  front  is  demonstrated  in  Fig.  11. 

The  pressure  and  particle  velocity  profiles  are  "rounded"  at  the 
shock  front  and  a  small  "tail"  precedes  the  actual  position  of  the  shock 
front.  This  typical  characteristic  of  the  numerical  method  produces 
an  inherent  "edge"  effect  in  the  calculations.  The  "tail"  causes  the 
unreacted  explosive  at  the  cylindrical  surface  to  expand  before  it  is 
"detonated"  by  the  computer  program.  Therefore  the  lateral  rarefaction 
wave  proceeds  into  the  cylinder  about  0.3  microseconds  sooner  than  it 
would  if  the  shock-detonation  front  were  infinitely  thin. 

The  specific  volume  (or  density,  since  p  =  1/V)  distribution 
in  the  cross  section  of  the  detonating  cylinder  at  t=11.3841  micro¬ 
seconds  is  shown  in  Fig.  12.  The  isochore  furthest  from  the 
detonation  front  is  V=1 . 50  cm3/g.  The  specific  volume  decreases  by 
0.05  cm3/g  on  each  consecutive  isochore  closer  to  the  detonation  front. 
At  the  detonation  front,  the  specific  volume  is  V=1/p  =1/2.210=0.4525 

cm3/g.  The  specific  volume  distribution  in  the  detonation  head  region 
is  identical  to  the  specific  volume  distribution  in  a  spherical  Taylor 
wave . 

VI.  DISCUSSION 

The  results  of  the  calculations  (Sec.  V)  are,  of  course,  only 
valid  for  that  part  of  the  detonation  where  our  physical  and 
mathematical  model  of  the  detonation  is  consistent  with  the  experi¬ 
mental  data.  More  specifically,  it  is  valid  over  that  region  where 
the  shape  of  the  detonatxon  front  can  be  represented  by  an  expanding 
sphere  of  radius  D  t.  From  the  extensive  experimental  data  of  Cook, 
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et  al1  , ,  this  is  shown  to  be  for  a  charge  length  of  3.5  charge 
diameters  for  pentolite,  in  particular,  and  for  charge  lengths  of  1 
to  4  charge  diameters  for  solid  explosives,  in  general.  Since  the 
calculations  were  carried  out  to  where  the  detonation  front  was 
1.71  charge  diameters  from  the  point  of  initiation,  we  are  well  with¬ 
in  the  limitations  of  our  model. 

In  addition  to  the  specific  numerical  results  obtained,  more 
general  conclusions  can  be  made.  We  have  shown  that  the  flow 
immediately  behind  the  detonation  front  in  the  region  unaffected 
by  the  rear  and  lateral  rarefaction  wave,  is,  identically,  the 
Taylor  similarity  flow  for  a  spherical  detonation.  Since  the 
characteristic  of  similarity  flows  is  that  the  flow  variables  depend 
only  on  their  relative  position  in  the  flow,  we  are  able  to  determine 
the  flow  in  this  region  for  later  times.  That  is,  at  least  to  the 
point  where  the  radius  of  curvature  of  the  detonation  front  is  equal 
to  3.5  charge  diameters  which  is  the  upper  limit  of  our  detonation 
model  according  to  Cook's  data.  From  the  rather  limited  experimental 
data  of  Jaffe  and  Clairmont4 ,  the  indication  is  that  our  model  of  the 
detonation  may  hold  over  a  much  greater  distance.  Therefore,  perhaps, 
the  Taylor  similarity  solution  may  be  validly  extended  to  charge 
lengths  greater  than  3.5  charge  diameters. 

It  should  be  pointed  out  that  although  one  can  obtain  the 
flow  in  the  detonation  head  region  for  later  times  because  it  is  a 
similarity  flow,  one  cannot,  however,  determine  the  shape  of  this 
region  by  extrapolating  the  calculated  results.  The  shape  of  the 
detonation  head  has  to  be  determined  by  calculating  the  motion  of  the 

lateral  rarefaction  front  into  the  spherical  Taylor  wave  flow  of  the 
detonation  head. 


13 


NOLTR  66-150 


VII .  ACKNOWLEDGEMENTS 

The  author  is  indebted  to  Dr.  H.  M.  Sternberg  for  suggesting 
the  problem  and  providing  helpful  advice  during  the  course  of  making 
the  numerical  calculations,  and  to  Drs.  D.  Price  and  S.  Jacobs  for 
taking  time  for  consultations  and  for  their  interest  in  the  problem. 
The  author  also  wishes  to  thank  H.  Hurwitz  for  making  the  RUBY 
thermochemical  calculations  and  W.  A.  Walker  for  making  available 
his  computer  routines  for  generating  the  Taylor  wave. 


NOLTR  66-150 


VIII.  REFERENCES 


1.  M.  A.  Cook,  et  al.,  J.  Appl.  Phys.,  21_,  269  (1956). 

2.  T.  Orlow,  D.  Piacesi,  and  H.  M.  Sternberg,  Naval  Ordnance 
Laboratory  Report  NOLTR  7265,  (September  1960) . 

3.  Ia.  B.  Zeldovich  and  A.  S.,  Kompaneets,  Theory  of  Detonation, 
(Academic  Press,  New  York  and  London,  1960) . 

4.  I.  Jaffe  and  A.  R.  Clairmont,  Naval  Ordnance  Laboratory  Report 
NOLTR  65-33,  (June  1965). 

5.  J.  Von  Neumann  and  R.  D.  Richtmyer,  J.  Appl.  Phys.,  2 1,  232  (1950). 

6.  J.  W.  Enig,  J.  Math,  and  Phys.,  Vol  XI,  No.  1,  (April  1961). 

7.  N.  L.  Coleburn  and  T.  P.  Liddiard,  Jr,,  J.  Chem.  Phys.,  44 
1929,  (1966). 

8.  H.  M.  Sternberg  and  D.  Piacesi,  Phys.  Fluids,  9,  1307,  (1966). 

9.  H.  B.  Levine  and  R.  E.  Sharpies,  University  of  California, 

Lawrence  Radiation  Laboratory  Report  UCRL-6815,  (March  1962) . 

10.  R.  D.  Cowan  and  W.  Fickett,  J.  Chem.  Phys.,  24,  932  (1956). 

11.  C.  L.  Mader ,  Los  Alamos  Scientific  Laboratory  Report  LA-2900, 

(July  1963). 

12.  G.  I.  Taylor,  Proc.  Roy.  Soc.,  A,  200,  235  (1950). 

13.  W.  A.  Walker  and  H.  M.  Sternberg,  presented  at  Fourth  Symposium  on 
Detonation,  (Oct  12-15,  1965),  to  be  published  in  the  proceedings 
(1967)  . 


15 


NOLTR  66-150 


FIG  .  1  A  QUALITATIVE  DESCRIPTION  OF  THE  FLOW  IN  A  POINT-INITIATED  CYLINDER  OF  EXPLOSIVE 


FIG.  2  A  GENERAL  LAGRANGIAN  COMPUTATION  GRID.  MASS  POINTS 
ARE  LOCATED  AT  THE  INTERSECTION  OF  THE  k  J  LINES. 
POSITIONS  AND  VELOCITIES  ARE  ASSOCIATED  WITH  EACH 
MASS  POINT.  PRESSURE,  SPECIFIC  VOLUME,  INTERNAL  ENERGY, 
AND  AN  ARTIFICIAL  VISCOSITY  ARE  LOCATED  AT  THE  CENTER 
(X's)  OF  THE  QUADRILATERALS  ("ZONES"). 


NOUR  66- 1 50 


OC  V 
UJ  TZ  “ 

uj  2 

fc 

8. 

Q  u. 


i  § 
-  x 

o  j8 

u- 

-c  P 


o 

_c 


CL 

O 


O  41 


< 

x 

u 


T> 

c 


>>• 

o 


<  u 

I-  E 

z  * 

“  *D 

■s.  _C 

^  u 

LU  .£ 

'i  % 

<  - 
t—  ID 

SH  u- 

o  ° 

i  I 

X  a 


N  U 
n 


v 

E 

E 


«.  o 

°  o 

8)  ii 

•i? 

W  N 

O  ' 

t  ^ 

o  e  e: 

t/i 

x  -c  o 
°  “  «, 

I  °  c 

0,?-° 
>  £  o 
Q  _o  X 
o  cl  ‘c 


•o  w 

S).£ 

2  8. 

*"  (U 

i*Si£ 


o 


(a313WVia  30HVHD  "1  VI1.I N i 


ijrivisia  iviova)  p/a 


NCLTR  66-150 


(8313WVIQ  39HVI-D  1VIJLINI  /  3^NV1SIQ  IVIQViO  P/H 


DETONATION -PRODUCT  GAS 


MOHR  66- 160 


/Hll.'WVICJ  lO'i 


•/ :  la  iviavai  p  a 


FIG.  5  Cross  section  of  the  detonating  cylinder  at  6.249  /Jseconds 
after  initiation.  The  detonation  front  has  reached  the  side  of 
the  cylinder  and  the  detonation-product  gas  expands  radially. 
Point  of  initiation  is  at  (R/d,Z/d)=  (0,0). 
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adjacent  isobar  lines  (actually  a  set  of  discrete  points)  with  the 
leftmost  isobar  =  1  Okb  and  the  pressure  at  the  detonation  front 
=  245kb.  Point  of  initiation  is  at  (R/d,Z/d)=  0,0. 
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The  pressure  distribution  in  the  form  of  isobcrs  is  shown  in  the 
cross  secrion  of  the  detonoting  cylinder  at  7.6163  fJ  seconds 
after  initiation.  There  is  a  difference  of  lOkb  between 
adjacent  isobar  lines.  Point  of  initiation  is  at  (R/d,Z/d)=  (0,0). 
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